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Abstract 

The problem of nuclear decays occurring due to the neutron-antineutron os- 
cillation is analyzed in the optical potential model and within the field-theoretic 
S-matrix approach. The result of the optical potential model for the nuclear 
decay width is rederived within the field-theoretic S-matrix approach. The n-ii 
oscillation in nuclei is suppressed drastically. We discuss an example for the 
spin precession of an atom, in which analogous suppression takes place due to 
transparent physical reasons. 
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The possible occurrence of the baryon-number-violating processes predicted by 
the grand unified theories (GUT's) stimulated searches for the neutron-antineutron 
oscillation 1-8 . The limits for the neutron oscillation time in the vacuum, r = 1/e > lyr, 
are derived from stability of nuclei 1-9 and experiments with free neutrons 10 . 

The n-ii oscillation is described by the Hamiltonian 

V = l -eJ dxtf (x)^* (x) (1) 
where o\ is the Pauli matrix and 

is a doublet of the neutron and antineutron fields. The V is Hermitian, since = 
( , n (x)) c = ?/> n (x) + (m the Majorana representation). 

For nuclei with atomic numbers A » 1, the decay width due to the n-ii oscillation 
takes the form 1-8 

r = e2jV (M-M) 2 + ry4' (3) 

Here N = A — Z is the number of neutrons in the nucleus, Z the number of protons, 
M the nuclear mass, and M the average mass of the state obtained from the nucleus 
by the replacement of a neutron by an antineutron. The value Ta is a hadron-scale 
antineutron annihilation width of the intermediate state. The nuclear decay rate per 
neutron, T/N = e£, contains an additional small parameter 

£ = < 10- 31 (4) 

* (M - M) 2 + ry4 1 ' 

as compared to the rate e of the vacuum n-ii oscillation. 

In recent paper 11 the finite-time S-matrix approach is developed for evaluation of the 
nuclear decay rate due to the neutron oscillation. It is claimed that the optical potential 
model is in disagreement with the field-theoretic S-matrix approach. We analyze more 
carefully the optical potential model, derive the Eq.(3) within the infinite- and finite- 
time S-matrix approaches, and clarify physical nature of the mechanism responsible for 
occurrence of the small parameter £. The potential model is shown to be in agreement 
with the S-matrix theory. 

We shall be interested in solutions to the evolution problem with the Hamiltonian 
of the form 
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H — ea± + U n P + + U n P- - iT{t)/2P_ (5) 

where P± — (1 ± 0" 3 )/2 are projection operators to the neutron and antineutron states, 
T(t) is a time dependent width. The real potentials, U n and U n , describe the neutron 
and antineutron rescattering effects on the surrounding nucleons. The potential model 
for the in-medium neutron oscillation is described by the Hamiltonian (5) with constant 
width T(t) = r^. The time dependent width T(t) is considered for further applications 
(see below). 

The Hamiltonian (5) can be split into two parts, H = H + V with H = U n P + + 
U n P- — iT(t)/2P_ and V = ea\. In the interaction representation, the potential V 
takes the form 

V I (t) = U- 1 (t t 0)VU(t,0) (6) 

with U(t,0) = exp(-i / H (t')dt') = e~ iUnt P+ + e^^'^Wp. and 7 (i) = \$T(t')dt'. 
The potential Vi(t) can be found to be 

Vj(t) = e(e- iAut -^a + + e +iAut+ ^a-) (7) 

where AU = U n — U n = M — M and a± = (<n ± i<r 2 )/2. 

The transition amplitudes are determined by the S-matrix elements in the Heisen- 
berg representation: S(t,0) = U(t, 0)Texp(— % JVj(t')dt'). To the second order in the 
e, we get 



S(t,0) = e- iU "\P + + e -* Aut -^P_ - ie{F_{t)a + + e ~ iAut -^ F+(t)a-) , . 

-e 2 (GL(t)P+ + e- lAut -<VG + (t)P_) + ...) [ } 

where 

G±(t) = Jo* e ±iAut,± ^F T (t)dt>. W 

The disappearance of the neutrons is described by a projection S(t,0)P + of the 
S-matrix to the neutron component of the initial wavefunction (2). Using Eqs.(8) and 
(9), we obtain for T(t) = T A and t » 1/T A 

S(t, 0)P + = e- iU ^(P + - lejxo- - eH/xP- + ...) (10) 
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where x = r a/2 + i(M — M). The decay probability is determined by absolute square 
of the S-matrix element | < i/j n \S(t, 0)\ip n > | 2 = 1 — e£t + where the use is made of 
the relation eRe{2/x} = £■ Comparison with the decay law for quasistationary states, 
| < Q\S(t,0)\Q > | 2 = exp(-r Q t) = l-T Q t + ... ait « 1/T Q allows to fix the nuclear 
decay width per neutron, T/N = e£. In this way we reproduce Eq.(3). 

The potential model operating with the non-Hermitian Hamiltonian (5) gives the 
phenomenological description for the nuclear decays. The results of the potential model 
can be justified within the field-theoretic S-matrix approach. 

Let 1 7V, Z> be a discrete energy eigenstate of the strong Hamiltonian Hq, associated 
to a nucleus of mass M, charge Z, atomic number A = N + Z, and zero momentum. 
With respect to the total Hamiltonian H = Hq + V (in the field-theoretic S-matrix 
approach, the Hq and V are Hermitian), the state \N, Z> is no longer an energy eigen- 
state. In what follows, the ceter-of-mass variables are factored out, so the normalization 
condition reads simply < N, Z\N, Z >= 1. The decay probability of the state \N, Z > 
to a channel X in the first order to the V is determined by the matrix element 

< X\V\N,Z >= eVNgxA(M x ). (11) 

The state \X > consists of N — 2 neutrons, Z protons, and a number of mesons. 
It belongs to the continuum energy spectrum of the Hq. The factor takes into 
account structure of the operator V. The function A(M X ) has the Lorentz form (cf. 
Ref.12, Ch.27) 

A 2 (M X ) = ^4, o-n- (12) 

v X! (M-M) 2 + ry4 v ; 

The state V\N, Z > has indefinite energy and definite (zero) momentum, since [H , V] ^ 
and [P, V] = 0. The value M is the average mass of this state. The scattering states 
\X > are grouped around the invariant mass Mx ~ M. The function A(Mx) reflects 
spread in the Mx of the scattering states X entering a superposition to form the lo- 
calized wave packet V\N, Z >. The state V\N,Z > decays with a strong width Ta- 
The value g x is a coupling constant to the channel X. This is a smooth function of 
the M x . 

In the nonrelativistic approximation, the state V\N, Z > consists of the fixed num- 
ber of particles. Its normalization condition reads < N, Z\V 2 \N, Z >= e 2 N . Inserting 
in this equation complete set of the states \X >, we obtain 



/ 



—A 2 (M')e(M') = 1 (13) 
2n 



4 



where 

e(M') = W g\{2-af5\P x - P')dr x , (14) 
x J 

P' = (M', 0), and dr x is element of the phase space. 
The expression for the nuclear width, 

r = E / I < X\V\N,Z > \ 2 (2n) 4 5\P x - P)dr x (15) 
x J 

where P = (M, 0) with the use of Eqs.(ll) and (14) takes the form T = e 2 NA 2 (M)e(M). 
For heavy nuclei T A « M and M - M « M. Substituting in Eq.(13) A 2 (M X ) w 
(2ir)5(M x — M) one obtains e(M) m e(M) 1. In agreement with Eq.(3), we get 
T ^ e 2 NA 2 (M). 

It is instructive to analyze the decay process in the time-dependent S-matrix ap- 
proach. In the interaction representation with respect to the strong Hamiltonian Hq, 
the state \N,Z > evolves according to the law \N,Z,t >= Si(t, 0)\N, Z, > where 
\N,Z,0 >= \N,Z > and 5/(4,0) = T exp(-i J Vj(t')dt'). The value V l is defined by 
Eq.(6) with U(t,0) = exp(—iH t). To the second order in the V, the decay amplitude 
takes the form 



< N, Z,0\N, Z,t>=l- C dh [ h dt 2 < N, Z\Ve~ i(H °- MKtl - t2) V\N, Z > . (16) 

Jo Jo 

Inserting complete set of the states X to the matrix element in the right hand side of 
this expression, we get 

< N, Z\Ve- l(H °- M ^V\N, Z >= t 2 N f —A 2 {M')e-^ M '- M ^e{M'). (17) 

J 2ir 

Here the use is made of Eqs.(ll) and (14). The values of the M' are grouped around 
the M by the Lorentz function A 2 (M'). Setting e(M') e(M) 1 and performing 
the integration, we obtain 

< N, Z\Ve-^ H °- M ^V\N, Z >^ e 2 Ne -(rA/2+i(M-M))t_ (lg) 
The decay amplitude (16) at t » 1/Ta takes the form 

< N,Z,0\N,Z,t>^ l-e 2 Nt/ X . (19) 
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The term 2Re{e 2 N/x} = tNl; appearing in the absolute square of this amplitude can 
be recognized as the nuclear decay width T. We derived therefore Eq.(3). 

The possibility of reproducing the potential model results within the field-theoretic 
S-matrix approach has recently been questioned 11 . We show that the results of Ref.ll 
are erroneous due to a misinterpretation of the optical theorem. 

The S-matrix determined by the Hamiltonian H = H Q + V can be expanded into 
a sum of the unit operator and the T-matrix: S — 1 — %T . The unitarity relation 
SS + = 1 implies 21 mT = -TT+ . 

The complete orthogonal basis \X >= {\X b >, \X SC >} can be constructed from 
eigenstates of the Hamiltonian Hq. The states \X b > are the finite-norm bound states 
(nuclei), while the states \X SC > are the infinite- norm scattering ones. The normaliza- 
tion conditions are < X b \X b >= 5bb>, < X SC \Y SC >= S(X — Y). The expansion of the 
unit takes the form 

1 = £|*6 >< Xb\ + / dfi(X sc )\X sc >< X sc \. (20) 
b J 

The value < X b \TT + \X b >with the help of Eq.(20) can be represented in the form 
< X b \TT + \X b >= W b + | < X b \(-iT)\X b > | 2 (21) 

where 

W h =J2\< X h \S\X b , > | 2 + / dfi(X sc )\ < X h \S\X sc > | 2 . (22) 
y^b J 

Here the orthogonality is taken into account: < X b \(—iT)\X b / >=< X b \S\X b > > for b' ^ 
b and < X b \(— iT)\X sc >=< X b \S\X sc >. The value W b is the total decay probability 
of the state X b . The second term in Eq.(21) can be transformed as | < X b \(— iT)\X b > 
| 2 = | < X b \S\X b > | 2 - 2Re < X b \S\X b > +1. The value | < X b \S\X b > | 2 = 1 - W b 
is the probability of finding the system t sec after its preparation in the initial state 
\X b >. Eq.(21) becomes 

< X h \TT + \X b >= 2(1 - Re < X h \S\X h >). (23) 

For the finite-norm states, the quantity < X b \TT + \X b > has no probability meaning, 
along with the imaginary part of the T-matrix element —21m < X b \T\X b >. 
The value < X SC |TT+|X SC > with the help of Eq.(20) is represented by 

< X SC |TT+|X SC >= W sc (24) 
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where 

W SC = J2\< X sc \S\Y b > | 2 + f dpL(Y sc )\ < X SC \S\Y SC > | 2 . (25) 
b J 

The term analogous to the second term in Eq.(21) does not occur because of zero 
measure of the infinite-norm scattering states. The value W sc is the total transition 
probability to the states Y ^ X sc . Eq.(24) is the optical theorem. 

The excited nuclei correspond to quasistationary states. These states represent 
localized wave packets and have finite norm like the bound states. One can show that 
for quasistationary states the relation (23) holds true with the substitution \Xb >^ 
\Q > ■ There is a deep analogy between properties of bound states and quasistationary 
states (see Ref.12, Ch.34). The parallel can also be drawn towards absence of the 
probability meaning for imaginary parts of the T-matrix diagonal amplitudes. The 
optical theorem takes place for the infinite-norm scattering states only. Notice that 
nuclei are bound states with respect to the strong Hamiltonian H and quasistationary 
ones with respect to the total Hamiltonian H including the baryon-number-violating 
part. 

It is worthwhile to make tree remarks, (i) The decay amplitude of a quasista- 
tionary state is given by the S-matrix element < Q\S\Q >. Up to a phase factor, 

< Q\S\Q >= e- r «*/ 2 , in which case < Q\TT + \Q >= 2(1 - e"^ 2 ). In Ref.10, the 
value < Q\TT + \Q > is misinterpreted as the decay probability and an expression 

< Q\TT + \Q >= 1 — e~ r °* is erroneously suggested. It is equivalent to the decay law 

< Q\S\Q >= (1 + e~ rQt )/2 which has nothing in common with the time behavior of the 
decay amplitude for quasistationary states. In particular, at t — > oo the quarter of the 
state \Q > survives. So, conclusions following Eq.(16) of Ref.ll are erroneous, (ii) It is 
claimed also that the infinite-time field-theoretic S-matrix approach is infrared diver- 
gent. In Fig.l of Ref.ll, the neutron line of the Green's function belongs to the nuclear 
Bethe-Salpeter wave function, and so the neutron is out of the mass shell. The decay 
amplitude therefore is not singular. Indeed, we showed that the infinite-time S-matrix 
approach (Eqs.(ll) to (15)) reproduces the result of the potential model (Eqs.(5) to 
(10)) and the result of the finite-time S-matrix approach (Eqs.(16) to (19)). (iii) The 
exponential decay law for quasistationary states takes place at 1/E < t < log(E/T)/T 
where E is energy release in the process (Ref.12, Ch.30). In our case E ~ 2 GeV 
and the value T is given by Eq.(3). The time interval for the exponential decay can be 
evaluated to be 1(T 25 < t < 10 38 sec for r = 1 yr, T A = 100 MeV, and N = 100. At an 
interval 1/E << t « 1/T, the decay probability is linear in time (e~ rQt = l—FQt+...). 
The quadratic time dependence 11 of the nuclear decay rate is in contradiction with the 
exponential decay law for quasistationary states. 
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To clarify physical nature of the decay rate suppression, we wish to discuss gedanken- 
experiment with deuterium ( 2 H) atoms in the spin j = 1/2 state, which exhibits 
transparent mechanisms for the spin-precession-rate suppression. 

Let us consider sequence of the Stern-Gerlach magnets (SG) which measure the 
spin- z projection and merge the outgoing atomic beams together (for a description 
such magnets see Ref.13). The SG magnets are placed along the x axis at equal 
distances Iq one from another, so that the atoms propagating along the straight line 
pass successively all SG magnets. The atoms entering the first magnet are polarized 
along the z axis. The existence of the magnetic field H = (H x , 0, 0) with a small 
component H x along the atomic beam (version I) is assumed outside the magnets. The 
atom spin precession is described by the Hamiltonian H = —ficrH with /i being the 
atomic magnetic moment (// ~ /is/3 where /is is the Bohr magneton). The atom has 
spin up at a time t with probability 

w T (t) =cos 2 (u; x t) (26) 

where uj x = —jj,H x . The neutron oscillation is described by the same equation with the 
substitution £ <-> u x . 

In the next experiment, the existence of the magnetic field H = (H x , 0, H z ) with 
H x « H z (version II) is assumed outside the magnets. The SG magnets are supplied 
also by detectors for the spin-down atoms. The time evolution of the spin-up atoms 
with the switched off detectors is given by equation 

WlI (t) = l-^sm 2 (cot) (27) 
ur 

where w = —fiH,uj = \uj\. If detectors are switched on, the evolution law (27) is no 
longer valid. The operating detector installed at the k-th SG magnet can produce no 
signal when a deuterium atom passes through it. In such a case, we know with cer- 
tainty that the atom has its spin up. The passage of the atom is accompanied by the 
wavefunction collapse. Afterwards, the atomic wavefunction up to a normalization fac- 
tor coincides with the initial one, and the oscillation process starts from the beginning. 
The probability that the first N detectors all give no signals is given by 

^//(t)on = (l-4 Sin2 M))' V ( 28 ) 

where N = t/t >> 1, to — lo/v is the time needed to pass from one magnet to the 
next, and v is the atom velocity. The time evolution of atoms in between the SG mag- 
nets is described quantum mechanically, whereas at t » t due to the measurement 
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procedure, the classical rules should be applied for calculation of the probability. Each 
instant of time t = kto the atom is either detected with a probability 1 — wu(to) or 
continue its propagation with a probability wji(t ). To derive Eq.(28), one should take 
N times the product of the elementary probabilities wu(t ). The value tujj(to) is close 
to unity, since uj x « u, so the probability to survive in the spin- up state follows the 
exponential decay law 

w n (t)on ~ exp(-u x £'t) (29) 

where uj x > is assumed and 

(30) 

to the first order in the uj x . 

In the first experiment, the value uj x determines precession rate of the spin- up 
atoms. In the second experiment, new dimensionless parameter £' << 1 occurs. The 
decay rate for the spin- up atoms turns out to be u x £' « uj x . 

The apparent analogy between the neutron oscillation and the atomic spin preces- 
sion is summarized below. 

neutrons < > spin-up atoms 

antineutrons < > spin-down atoms 

vacuum n-ii oscillation < — > atom spin precession (v. I) 

vacuum oscillation rate e of neutrons < — > spin precession rate u x of atoms (v. I) 

AU — M — M < — »■ u z (v. II) 

ii length free path in nuclear matter < — > distance l between SG magnets (v. II) 

antineutron annihilation < > detecting spin-down atoms (v. II) 

suppression parameter ^ < > suppression parameter (v. II) 

nuclear decay rate per neutron < > atom spin- up decay rate uj x ^' (v. II) 

Reasons, according to which physical results in these two problems are identical, can 
be explained by considering the atomic spin dynamics described by the Hamiltonian 
(5) with the time dependent width 

T(t) = J2^(t-kt ). (31) 

k 

Here the k's are integer numbers. In the intervals kt < t < (k + l)to the atom 
is out of the SG magnets and its spin precession is governed by the magnetic field 
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H = (H x , 0, H z ). At the time t = kto the act of the measurement takes place and 
the spin-down component of the wavefunction vanishes. In order to set the lower 
component equal to zero, we must pass to the limit a — > oo. The upper component 
at t — kt remains continuous. The probability to find the deuterium atom in the 
spin-up state can be computed according to the usual rule as the modulus squared 
of the upper component of the atomic wavefunction. Using such a prescription, the 
results of Eqs.(28) and (29) can easily be reproduced. 

One can show that the S-matrix projection S(t,0)P+ to the spin-up initial state, 
being averaged over an interval (t — T,t + T) at T >> to, 



coincides with the S-matrix projection for the constant width Ta- It means that projec- 
tion S(t, 0)P + responsible for disappearance of atoms from the beam can be reproduced 
in a model identical to the potential model for description of the n-n oscillation in nu- 
clei. 

In the limit a — > oo, the functions entering Eq.(8) can easily be computed. After 
averaging over the time, we get the following expression 



where 1/x' = {sin 2 {uj z to) — i((uj z to) — sin{ui z to)cos{ui z to))) / {2oj 2 to) ■ Notice that u x Re 
{2/x'} = Eq.(33) is identical to Eq.(10) for e — u x and x — x' ■ Given that Ta and 
AU are known, one can find to and uj z . To the second order in the uj x (e), the evolution 
laws in these two problems are identical in the average sense. 

The qualitative explanation for the neutron-decay-rate suppression can be given in 
the following way. In the case AU >> Ta one can start from Eq.(31) that gives average 
admixture ~ w 2 /w 2 of the spin-down atoms in the beam. The rate of disappearance 
of the atoms can be estimated to be ~ u) 2 /{uj 2 to) « uj x . The equation x — x' 
gives uj z = AU/2, t ~ l/T A >> 2/ AC/. Using the above analogy, we derive the 
corresponding neutron decay rate e 2 T A /AU 2 = << e. 

The inequality AU ^ is not unique reason for suppression of the decay rate. Let 
AU = 0. The equation x — x' gives uj z = and to = 4/T^. The decay rate of the 
spin- up atoms becomes u x ^' = uJ x t « uo x . The detection of the spin-down atoms 
suppresses the spin precession process by itself. The neutron propagation in nuclei is 
accompanied by interactions with surrounding nucleons. The meson exchanges take 
place each t = 4/T^ 1CT 24 sec. They can be interpreted as acts of measurements 




(32) 



S(t, 0)P + = e^\P + - icujx'a- - ui 2 x t/x'P + + ...) 



(33) 
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of the lower antineutron component of the wavefunction (2). These interactions are 
followed by collapse of the wavefunction to the pure neutron state. The wavefunction 
collapse in turn reduces the nuclear decay rate down to Ae 2 /T A = e£ << e. 

The analogous suppression mechanism exists for ultracold neutrons in a trap. The 
absorption probability for antineutrons in collisions with the copper boundary is estimat- 
ed 3 to be ~ 1/5 for the tangent velocities ~ 4m/ sec. The suppression of the precession 
rate occurs due to the wavefunction collapse that takes place with certainty under those 
conditions one time per 5 neutron collisions with the trap boundary. The suppres- 
sion factor is about 5et < 10~ 6 , where to (~ seconds) is a typical time between two 
collisions. 

In the conclusion, the field-theoretic S-matrix approach allows to justify the results 
of the potential model for heavy nuclei. The suppression of the nuclear decay rate 
is the result of the energy gap AU and the antineutron annihilation width Ta- The 
suppression of the neutron oscillation process is not specific for nuclear physics. In the 
SG-type experiment described above, the atom spin precession is suppressed due to the 
similar physical reasons. For light nuclei, the procedure 4 ' 5 of solving the inhomogeneous 
Schrodinger equation for antineutron wavefunction in the optical potential is more 
adequate then the direct use of Eq.(3). Being applied for heavy nuclei, such a procedure 
reproduces Eq.(3). The models 1-7 account to a various degree of accuracy for dynamics 
of the antineutron creation and annihilation in nuclei. Further refinement of parameters 
of the antineutron optical potential entering Eq.(3) would be desirable. It is unlikely, 
however, that the corresponding numerical estimates are changed more then by an 
order of magnitude. 

After submitting this work our attention has been called to a recent paper [14] in 
which results of Ref.ll obtained with the use of the field-theoretic methods are critisized 
on the basis the optical potential model. The authors give two new derivations for 
expression for the nuclear decay width within the optical potential model (Eqs.(l) to 
(5) and Eqs.(ll) to (14)). We argued that the imaginary part of the T-matrix has no 
probability meaning for bound and quasistationary states. It is not clear then why 
the value W n in Eq.(3) of Ref.14 is not a probability whereas the value W in Eq.(5) 
does. These values both are defined, respectively, as imaginary parts of strong and 
baryon-number violating (GUT's) T-matix amplitudes. We discussed the T-matrix 
properties for Hermitian Hamiltonians. The optical model Hamiltonian H = —iT/2 
used in Ref.14 is not Hermitian, and so it is necessary to reanalyze the problem on 
slightly more general grounds. 

Notice that Eq.(16) of Ref.ll which represents afinite-time analog of Eq.(4) of Ref.3 
is correct and it remains valid for non-Hermitian Hamiltonians also, since hermiticity 
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of the Hamiltonian is not assumed in deriving this equation. It is allowed therefore 
to substitute here H = —iV/2. In this way, purely imaginary T-matrix occurs: T = 
-iW/2 with W = effort » 1/T A . The S-matrix becomes S = 1-iT = l-W/2. The 
probability of finding the system in the initial state t sec after its preparation becomes 
\S\ 2 = 1 + 2Re(-iT) +TT+ w 1 + 2ImT = 1-W. The quantity W = -2ImT thus 
really receives, to the first order in T — 0(t), a meaning of the decay probability. We 
can expand the value W n in Eq. (3) of Ref.14 in power series in t(— t a — tp) and verify, in 
agreement with that conclusion, that the lowest order term gives the decay probability 
Fa t « 1. The reason stems from the fact of neglecting the small term TT + = 0(t 2 ) 
in the expression for the IS*) 2 . At high times, however, W n — > 2. This is no longer 
surprise, since TT + becomes high as t increases, and so it should be taken into account 
in the IS*) 2 . The value W defined by Eq.(l) in Ref.14 is calculated to the first order in t. 
This is a reason for it gives the correct result for the probability. The next order term 
in t, however, should already be wrong, if the value TT + is disregarded. The instability 
of quasistationary states implies < Q\S\Q >^ when t increases. Our Eq.(23) gives 
then W n -> 2 at t » l/T A and W -> 2 at t » 1/T. 

Therefore, the value W n in Eq.(3) of Ref.14 is not a probability, when it is computed 
to all orders in t. To the first order in t it acquires the probability meaning. The value 
W in Eq. (5) is not a probability also. It coincides with the probability to the first order 
in t only. 

The author is indebted to K.G.Boreskov, B.V.Martemyanov and M.G.Schepkin 
for useful discussions, A. Gal for communication on the paper [14], and especially to 
L.A.Kondratyuk for attracting the interest to the problem of the neutron oscillations 
and valuable discussions of the optical potential model. This work is supported by the 
Russian Fund for Fundamental Researches under the Grant No. 94-02-03068. 
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